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On In-Plane Current Distribution Producing a Given Axisymmetric 
Distribution of Normal-to-Plane Magnetic Field 

Yu. E. Kuzovle^Q 

Donetsk Physics and Technology Institute, 83114 Donetsk, Ukraine 

A very good as well as simple approximate solution is suggested to the problem which flat circular 
distribution of electric current in finite film produces a given axially symmetric distribution of 
normal-to-film magnetic field. 
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1. In films of type II superconductors, a concentration, 
n , of vortices piercing the film usually simply correlates 
with magnetic flux density, B , at its surface: ^oi^ — 
(B — Bo)n , where ()„ means normal-to-film component 
of a vector, $o is magnetic flux quantum, Bq - uniform 
external field, and B — Bq is magnetic field induced 
by (super-)current flowing in the film, J . The same 
current stimulates the vortices to move. To consider this 
motion, analytically or numerically, we should be able to 
determine what is distribution of J which produces a 
given distribution oi B — Bq = ^qu . 

Such task is non-trivial even under flat geometry, if 
the fllm physically occupies a part of its plane and cor- 
respondingly the field B is considered to be known in 
that part only. Previously, an exact analytical solution 
of the task was found for the single case of infinitely long 
flat strip [3,13, assuming homogeneity of B along it. 

Recently, I suggested exact solution for one more case 
of round fllm j3| , assuming axial symmetry of B . The 
obtained radial distribution of circular current, J(r) , as 
functional of radial distribution of H{r) = {B — So)n , 
can be represented by sum of regular term and singular 
term: 



Jir) = JREcir) + JsiNcir) 



jREG{r) = -( T{r, p)p '^^j-P^ dp 
Jo 



dp 



JsiNcir) = 



' H{p)pdp 



(1) 



(2) 



(3) 



(4) 



Here, J has the sense of the current density, j , multi- 
plied by 27r/c (in CGS units) and integrated over fllm's 
thickness; fllm's radius is taken to be unit of length; the 
kernel T{r,p) is deflned by formulas 

TT min(r, p) 



k = 



min(r, p) 
max(r, p) 



= arcsin max (r, p) , 



where standard designations of the complete and incom- 
plete elliptic integrals of 1-st and 2-nd kinds are used (see 
e.g. U]). I permit myself to copy from [3] the plots of 
T{r,p)p (see Fig.l below). 

Notice, first, that the singular contribution JsiNciT) 
looks exactly as Meissner current in vortex-free fllm un- 
der effective uniform external fleld H^, . Second, the reg- 
ular part of the current, JREcif) , turns into zero (or 
at least remains finite) at fllm's edge. Thus the condi- 
tion He ~ ensures the absence of any current singu- 
larity at the edge. At that, the relation between J(r) 
and H{r) becomes non-local analogue of the local one, 
Anj /c — ~dH/dr , which takes place under cylindric ge- 
ometry. 

2. In general, it may be convenient to unify Jreg{t) 
and Jsing{i") into the single expression Q: 



J(r) 



G{t,p) 



k = 



d 
dr 



n {p + r) 



G(r,p)pH(p) dp, 



[K{k)~2Fi<t>,k)] , 



(6) 



(7) 



pr 



= — ( arcsin p + arcsin r ) 



Both © and (0) contain incomplete elliptic integrals. 
But, unfortunately, not all of popular mathematical pro- 
grams can quickly (or let somehow) calculate them. This 
fact just implies the subject of present article, namely, 
approximation of the kernel G(r, p) in terms of complete 
elliptic integrals only. 

For this purpose, we can apply the series expansion of 
the kernel G{r,p) (see 0): 



G(r, p)^'^Y.{Ak + 3) Si P2k+i (/) P2k+i [p') (8) 

TT ^ — ' 

fc=0 



where x' = vT 
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Sk = B(l/2,fc+l)/2 = 2'=/c!/(2fc-f 1)!! , (9) 
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FIG. 1: T{x,y)y via y at x = 0.1, 0.3, 0.5, 0.7 and 0.9. 




which can be easy proved with the help of the Stirhng 
formula. 

Our approximation will be merely replacement of the 
factor (4fc + 3)5| by its lower bound, that is by 3 , in ac- 
cordance with H1U|I and © • Corresponding approximate 
series for the kernel G{r, p) is 



G{r,p) « Gapp{r,p) ^ - J2 ^2fe+i [r') P2k+i ip') 

fc=0 

(11) 

Fortunately, its sum can be fully expressed via complete 
elliptic integrals, if use the representation 

oo ^ 

P2k+i{x)P2k+i{y) = -[F{x,y)~F{x,~y)] , 



k=0 



F{x,y) 



where P(<, u) is generating function of the Legendre 
polynomials. 



p{t,u) ^Y.ep^{u) 



1 



n=0 



Vl + <^ - 2tw 



and then reduce F(x^ y) to standard elliptic-type inte- 
grals HH]. The result is 



G app {•^i y) 



3^/2 



Kik^) K{k+) 



(12) 



FIG. 2: G(x, y) y (pink curves) in comparison with 
Gapp{x,y) y (black curves) via y at x = 0.1, 0.3, 0.5, 0.7 
and 0.9 . 



with B{x,y) = r{x)T{y)/r{x + y) being the beta- 
function, and P„i are ordinary Legendre polynomials Q . 
The key to the desired simplification is in the inequalities 



1 < 



4k + 3 
12 



1.0472 , (10) 



S±^l + xy±^il-x^)il-y^) , fc±^^ (13) 

Fig. 2 shows this approximation in comparison with for- 
mally exact kernel (in your mind please continue the log- 
arithnic peaks to infinity). In practice, relative error of 
the approximation does not exceed 3 % . 

I am grateful to Dr. Yu. Genenko and Dr. S. Yampolskii 
for useful comments and attracting my attention to the 
work 0. 
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